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In this paper we develop a new approximation method 
valid for a wide family of nonlinear wave equations of Nonlin- 
ear Schrodinger type. The result is a reduced set of ordinary 
differential equations for a finite set of parameters measur- 
ing global properties of the solutions, named momenta. We 
prove that these equations provide exact results in some rele- 
vant cases and show how to impose reasonable approximations 
that can be regarded as a perturbative approach and as an 
extension of the time dependent variational method. 

PACS: 03.40.Kf, 42.65.-k, 03.75.-b 



The solution of nonlinear wave equations representing 
physically relevant phenomena is a task of the highest in- 
terest. However, it is not possible to find exact solutions 
except for a few simple cases where one is lucky to inte- 
grate the equations involved. In particular, in the 70's 
some mathematical techniques were discovered which al- 
lowed the integration of several relevant nonlinear wave 
equations ffl. So the development of rigorous approxi- 
mation methods is of interest in those cases where the 
equations are known to be non-integrable. 

One family of nonlinear wave equations with lots of 
practical applications is that of Nonlinear Schrodinger 
Equations (NSE) [pH, which arise in plasma physics, 
biomolecule dynamics, fundamentals of quantum me- 
chanics, beam physics, etc., but specially in the fields 
of Nonlinear Optics j| and Bose-Einstein condensation 
H . In the last two fields a great variety of these equa- 
tions appear involving different spatial dimensionalities, 
nonlinear terms (saturation, polynomial, nonlocal, losses, 
etc.) and number of coupled equations. 

One common approximate theoretical approach to the 
analysis of the dynamics involved in those problems is 
to assume a fixed shape for the solution with a finite set 
of time dependent parameters so that the dimensional- 
ity is made finite at the cost of loss of information on 
the solution (many degrees of freedom are lost). This 
method receives many denominations depending on the 
context: collective coordinate technique, time-dependent 
variational method, equivalent particle approach, energy 
balance equations, etc. 

Although not explicitly stated most of those methods 
can be reduced to a more elegant formulation which is 
the time dependent variational technique, originally de- 
veloped by Anderson |(| for one dimensional problems 
based on Ritz's optimization procedure. This approxi- 



mate technique is a good tool to study the propagation 
of distributions having simple shape. If the shape of the 
actual solution is close to the trial function, the outcome 
of variational method will be in good agreement with the 
real solutions, otherwise it may be very rough or even fail 
0. Despite of this fact the technique has been used in 
many physical situations. In Nonlinear Optics it has been 
applied to many problems some of them being listed in 
Rcf. The method has been applied to many other 
physical problems where nonlinear wave equations (in 
particular NSEs) arise including random perturbations 
H , nonlocal equations |[o| , collapse phenomena jO,^2[ . 
propagation and scattering of nonlinear waves |13||, etc. 
A review of the application of the technique with em- 
phasis on problems with different scales (focused in con- 
densed matter) is given in Ref. Q]. This technique has 
been also used in the last years in the framework of Bose- 
Einstein condensation (BEC) applications to explain the 
low energy excitation spectrum of single |ll| and double 
Jl2| | condensates, collapse dynamics [|l2|, and many other 
problems ]l6||r^ ], etc. 

In this letter we develop a completely different tech- 
nique called the moment method [Q. This method is 
based on the definition of several integral parameters 
whose evolution can be computed in closed form and has 
been used to obtain exact results in particular applica- 
tions p9|-pl|. We provide here a general framework for 
its application as well as several ways to treat it system- 
atically as a perturbative technique. 

The general NSE and moment equations.- Let us con- 
sider the rt-dimensional NSE 



V(r)^ + g(\^\ 2 ,t)^-ia(\^\ 2 ,t)ij (1) 



where A = £ fc d 2 /dx 2 , and V{r) = J2 k ^ k {t) 2 {t)x 2 is a 
time and spatially dependent parabolic potential which 
has been included because it is always present in BEC 
problems. To aid readability, we will separate the so- 
lution in modulus and phase, ip = ^fpe 1 ^ , and define 
an interaction energy density G(p) as g = dG/dp. The 
nonlinear terms will be analytic functions of p such that 
g(p),G(p)->Q, p^O. 

In Table | we define the so called momenta of ip. Some 
of them are related to the momenta of the distribution 
p = \i/j\ 2 , and they all have a physical meaning (see Ref. 
p9[ for their interpretation in Optics). The evolution 
equations for these quantities are 
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Here DG is a shorthand notation for G(p) — g(p)p. 
Through this paper we will concentrate in the most com- 
mon case a(\ip\ 2 ,t) — a(t) for which the equations are 
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As stated before, some of these laws can be found in other 
treatments p|]l9[] which mostly concentrate on particular 
cases or treat them as basis for perturbation methods 
choosing one particular shape for the solution. From Eqs. 
(pT]-p(]) we find exact closed equations for zeroth and first 
order momenta, 



d 2 X l 
dt 2 



-u i (t)X i -2a(t)^-2a(t)X i . (4) 

In conservative systems with any type of potential 
the classical result of Quantum Mechanics d(xi) /dt — 
— (dV/dxi) is obtained as described in JlT[ | - 

Once Eqs. ( |3a||3"c| ) are integrated one is left with the 
problem of solving the remaining 3n + 1 equations. Typi- 
cally these equations do not form a closed set but involve 
integral quantities which are not included in the defini- 
tions of the momenta. To close them, i.e. to equal the 



numbers of equations and of unknowns, one must either 
restrict the problem or impose some kind of approxima- 
tion. 

Exact closure of moment equations.- We have found 
only two relevant cases in which the closure is exact for 
the rest of the equations. Both simplified problems cor- 
respond to conservative, a — 0, spherically symmetric 
potentials, Ui = uj(t) where 
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being D — J[G — gp] and R = VW the radial width. The 
first integrable case corresponds to n — 2 and G = Up 2 . 
For it Eqs. (ph simplify to 



d 2 R , s M 

- w = -.{t)R+ W 



(6) 



Here M a constant that depends only on the initial data 
and interaction strength U . This equation has been used 
to prove the existence of extended resonances in Ref. . 

Another ample family of systems for which moment 
equations are closed and exact is those with a time- 
independent interaction strength, dG/dt = (the usual 
case) , and a divergenceless velocity distribution (given by 
the phase gradient) div (V0) = A<fr = 0. This condition 
imposes no restriction on the density distribution p and 
is automatically satisfied by the well known vortex-line 
solutions, which in n spatial dimensions read 



ij> = p B (xi, ...,x n ;t)e l<l,B , 

Xk - Xk 



t>B = aX + arctan 
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Here X are free, time dependent parameters, and po is 
arbitrary. With those simple conditions we get an in- 
finite number of constants of evolution named "super- 
momenta", Q{F) = J F(p), built up from differentiable 
functions of the modulus, F(p), that satisfy the regular- 
ity conditions. In these cases one can prove that D is a 
constant and that Eqs. (||) become equivalent to Eq. (|^). 

Uniform divergence approximation.- Intuition dictates 
that the zero-Laplacian phase condition is related to (a) 
the configuration of the cloud, i.e. p, does not change, 
(b) the soliton or the wavepacket is either stationary 
or at most suffers displacements and rotations and (c) 
all of the "supermomenta" depend solely on the norm, 
Q(F) — f(N). General solutions have nonzero diver- 
gence of the velocity field, thus the zeroth order approx- 
imation, Acj) — 0, fails to describe the system. The 
next possibility is a first order approximation in which 
the Laplacian of the phase is uniform. As we will see, 
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this extension now allows for changes in the shape of the 
cloud and introduces three new independent variables in 
the supermomenta Q(F). Mathematically, the first order 
approximation the phase is 



4> = 4> B (x, t) + ya J (t)x 



(8) 



where <pg(x,t) is any function satisfying = 0. This 
approximation will be called uniform divergence approx- 
imation on the phase in what follows. A limited ver- 
sion of this approximation that uses a linear function 
(which has zero divergence) in place of <pB was first ap- 
plied to radially symmetric problems in Ref. pC| ] and 
to the study of resonances in general 3D forced NSE 
problems in Refs. |]l7| , ^l| . In that case, were 4>(x,t) — 
<po{t) + ax + ^2j (3j(t)Xj it is evident that the approxima- 
tion consists on a Taylor expansion of the phase or even 
better a polynomial fitting with time dependent param- 
eters. Here we provide a general framework for the ap- 
plication of the technique to arbitrary NSE problems. In 
a certain sense this is a generalization of the usual time 
dependent variational method but now no assumption 
on the shape of the amplitude of the wave is needed and 
the phase has a free is approximated by a least-squares 
type fitting with time-dependent parameters and only 
very general restrictions on the phase. 

Let us first take the case with <j)B = ax. We will 
assume for simplicity that the nonlinearity is or can be 
approximated by a polynomial 



G( P ) = ]Ta fc p fc 



(9) 



The dissipative terms can be removed from the equations 
by rescaling the solution with 7 = J cr(t')dt' , so that 
p — e 7 p. We will denote the momenta obtained using p 
by a tilde, as in Wi. It is important to stress that Eqs. 
(||) do not close immediately and it is necessary to study 
the monomial supermomenta 



Their evolution laws are 
dQ^ 
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Since the parameters in the phase can be expressed as 



A = -iogV^ 



(12) 



one obtains a closed form for each of the supermomenta 

-ro+l 



g(m) = C (m) (^Jw 1 ...W d ^ 



(13) 



where the constants must be determined from ini- 

tial data. We can use this property to estimate all the 
integrals in which G{p) appears, as a function of pow- 
ers of the mean square widths. Using this, and defining 
as before the natural widths, R — \fwi, in the general 
nonsymmetric case we arrive to 

j2 6 

= -UiRi + — H - (14) 



Mi 



dt ' • Rf i?! 

where Mi is a constant to be determined from the initial 
data, and Q(R\ . . . Rd, t) is a function that may be calcu- 
lated from Eqs. (||) and ([T3|). Finally, to interpret these 
results one must remember that the actual width of the 
cloud is actually Ri(t) = R i (t)e~ 2 ~ f . 

Eq. ( |l4| ) means that given an initial data it is possible 
to compute the evolution of the width (and all the mo- 
menta of the initial datum) provided the parabolic phase 
is a good description for the solution. The method is very 
powerful as it accounts for the evolution of any solution 
that can be described in terms of a finite number of mo- 
menta (since higher order momenta are functions of lower 
order ones). Though much more general than the collec- 
tive coordinate method the one presented here is simpler 
to apply and extend since it only involves computing the 
integrals in Eqs. (||). 

We must also remark that the phase is not restricted 
to be a polynomial in r as in Ref. [^0| . Instead one still 
obtains information when the solution has many other 
forms. For example, if phase is of type (0), then one 
can combine the two widths from the plane on which the 
vorticity is present, Wk and Wi, into a radial one, R = 
Wk + Wi and the evolution of R can be studied together 
with the widths from the remaining spatial directions. 

Independent moment approximation.- We have seen 
that the uniform divergence condition leads to a closed 
expression for every Q( m ) in terms of the widths Wi. Al- 
though this approach is more powerful than the usual 
time dependent ansatz it is possible to improve the mo- 
ment method to higher precision. The basic idea is to 
assume that only a finite set of momenta are indepen- 
dent, the number of those independent momenta being 
related to the accuracy of the solution and finding ex- 
pressions for higher order momenta as functions of lower 
order ones. For instance if only the {N, Wi} momenta 
are truly independent (as it occurs with the uniform di- 
vergence approximation) it is possible to find expressions 
for the rest of the magnitudes in terms of these momenta. 
By scaling the wave with respect to one of the coordinates 

1 x- 
i/j(x 1 ,...,x i ,...) -> -==V>(xi, . . . , — — , . . .), (15) 



l + e 



and relating the first order changes in Q( m ) and Wk one 
arrives to 
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This equation has a solution of the form reflected in Eq. 
( |l3| ) and also similar expressions may be derived for the 
rest of the unknown integrals in Eq. (^a|). Thus, in the 
first order case discussed here the method can provide a 
closed set of equations for any type of nonlinear term. 
In principle the procedure could be extended to higher 
precision approximation of the evolution of initial data. 
Work on this point is in progress and will be reported in 
detail elsewhere. 

Finally we would like to point out that Eqs. (||) can 
be used straightforwardly by replacing ip{r) with an ap- 
propriate ansatz. In the conservative case, a = 0, this 
procedure is equivalent to Ritz's optimization procedure, 
with the advantages that one needs not build a huge La- 
grangian integral, and that our central equations (|2a| ) 
can also be applied to dissipative systems that lack a La- 
grangian density at all. It must also be remarked that 
unless the ansatz has a phase different from (||), one will 
always arrive to Eq. (pl|). 

In conclusion we have presented the moment method 
in a general framework and discussed under which condi- 
tions it leads to closed equations for a finite set of parame- 
ters. The uniform divergence ansatz has been introduced 
as a way to improve usual collective coordinate methods 
and obtain a closed set of equations for any NSE with 
polynomial nonlinearity and linear dissipation (with or 
without a parabolic time dependent spatial potential). 
The second method takes the zcroth, first and second or- 
der momenta to be the only independent ones and builds 
analytical approximations for other magnitudes. It is an 
extension of the uniform divergence approximation that 
seems to work for non polynomic self-interaction energies. 

This work has been partially supported by CICYT un- 
der grant PB96-0534. 
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TABLE I. Definitions of the momenta and physical inter- 



pretation. 


Definition 


Interpretation 


N = fp 

Wi = J xip 

J = JG( P ) 


Norm, number of particles or intensity 
Center of mass or of energy 
Speed of center of mass 
Widths 

Speed of growth 

Kinetic energy 

Self interaction energy 
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